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We analyse the mechanism of quintessential inflation in brane world scenarios for a number of par- 
ticle physics inspired scalar potentials. We constrain the parameter space of those scalar potentials 
and comment on the likelihood that we could discriminate these models from standard inflation, 
based on upcoming large scale structure and cosmic microwave background observations. 
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I. INTRODUCTION 

There has been interest, recently, on the cosmologi- 
cal implications arising in a certain class of brane world 
scenarios where the Friedmann equation is modified at 
very high energies [1]. The key feature is that the mod- 
ification makes it easier to obtain inflation in the early 
universe, by contributing extra friction to the scalar field 
equation of motion [2-4]. Analysing this situation, it 
has been explicitly demonstrated how inflation naturally 
occurs in such models even when the potential is too 
steep to allow conventional inflation to occur [5,6]. In- 
flation ends when the brane world corrections begin to 
lose their dominance, and reheating takes place as a re- 
sult of gravitational particle production, rather than the 
usual inflaton decay mechanism. In this case the inflaton 
energy density subsequently redshifts sufficiently quickly 
that the produced radiation comes to dominate. One of 
the nice features this allows for is that the inflaton po- 
tential does not have to have a feature in it, and the field 
does not have to decay completely as we do not require a 
minimum to the potential about which the inflaton nor- 
mally oscillates as it decays. Instead, we can consider the 
intriguing possibility that the inflaton actually survives 
as inflation ends. As the universe evolves through radi- 
ation to matter domination the inflaton potential does 
not play much of a role, however, at a recent redshift of 
z ~ 1, it once again comes to dominate driving the uni- 
verse to an accelerated expansion. This mechanism has 
been named "quintessential inflation" scenario [7,8]. In 
this paper we investigate this possibility in more depth. 
Similar approaches have been adopted in the brane world 
context [9-11] although in this paper we will be analysing 
a broader class of potentials and discussing the likelihood 
of being able to detect the presence of such a field in the 
next generation of cosmic microwave background (CMB) 
observations. The class of potentials we will investigate 
include inverse power law, exponential and potentials mo- 
tivated by supergravity models. 

We consider the five-dimensional brane scenario, in 
which the Friedmann equation is modified from its usual 
form, becoming [1] 
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where k 2 = 8irG — 8tt/M 2 and M 4 is the four- 
dimensional Planck mass. A4 is the four-dimensional 
cosmological constant and the final term represents the 
influence of bulk gravitons on the brane. 

The brane tension Ab relates the four and five- 
dimensional Planck masses via 
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and is constrained by the requirement of successful nu- 
cleosynthesis as Ab > (lMeV) 4 . We assume that the 
four-dimensional cosmological constant cancels out by 
some (as yet undiscovered) mechanism, and once infla- 
tion begins the final term will rapidly become unimpor- 
tant, leaving us with 
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We assume that the scalar field is confined to the brane, 
so that its field equation has the standard form 
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Following Ref. [12] the density of particles produced 
after inflation is 
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(5) 



where g p is the number of scalar fields involved in particle 
production, likely to be 10 < g p < 100, and Hend is the 
value of the Hubble parameter at the end of inflation. 

In our brane world scenario, the ratio of the energy 
densities of the scalar field to the relativistic particles is 
well determined from Eq. (5), 
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Hence, the ratio of the energy densities is inversely pro- 
portional to the number of fields involved in gravitational 
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production g p . Moreover, it decreases as the time be- 
tween the end of inflation a en d and the maximum of the 
energy density of the produced relativistic particles a p 
increases. Defining A = log(a p /a on d), we find from our 
numerical simulations that this value is weakly depen- 
dent on the parameters of the scalar potential or the 
number of e-foldings of inflation, and spans the interval 
0.29 < A < 0.4. Therefore, when constraining the pa- 
rameter space of scalar potentials, we will consider both 
the average case of A = 0.35 with g p = 100, and the 
case of instant reheating A = with g p = 10, which 
corresponds to the situation when the quoted ratio is its 
maximum. 



II. INFLATIONARY PARAMETERS 

It proves useful to define the slow-roll parameters, 
analogously to the usual inflationary case. Following 
Refs. [3,9] we define e, r\ and £ 2 , generalising the usual 
ones [13], by 
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where prime indicates a </>-derivative and the slow-roll 
approximation has been employed. 

As the Universe inflates, the associated number of e- 
foldings of the scale factor a, is given by 
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where a cn d corresponds to the end of inflation, given by 
the condition e = 1. In this model, the number of e- 
foldings corresponding to scales entering the Hubble ra- 
dius today can be unambiguously determined [11]. The 
scale that leaves the Hubble radius during inflation and 
re-enters today is k = clnHn = aoH , i.e. 
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The subscript N corresponds to the value of the quan- 
tity N e-foldings from the end of inflation. Noting that 
aw/aend = cxp(-TV) and p^ nd = p^ (ao/«p) 4 we find 
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where = 2Al\h~ 2 x 10~ 5 is the fractional energy 
density in radiation today. 

The amplitude of the scalar and tensor perturbations 
produced during inflation are given by [3,9] 
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In the low energy limit, p <C A b (x <C 1), F 2 
whereas F 2 — > 31 / /2A b in the high energy limit. 

Having obtained the amplitude of the scalar and tensor 
perturbations the corresponding scalar and tensor spec- 
tral indices can be obtained. Using 
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we find [9] 
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If we define the quantity S to be for the standard 
inflationary scenario (without the brane corrections) and 
5 = 1 in the brane inspired case, then we can simplify 
Eq. (19) for the low and high energy limits as, 
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We also obtain that the consistency equation is the same 
in either limit [9], 
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The derivatives of the scalar and tensor spectral indices 
give in the limits of high and low energies 
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In the next section we turn our attention to a num- 
ber of particle physics inspired potentials, to determine 
how the brane corrections modify the dynamics of the as- 
sociated scalar field, leading to inflation in a number of 
scenarios where inflation would not have normally been 
present. In the following section we will proceed to in- 
vestigate the class of models where the inflaton field then 
becomes the quintessence field by today. 
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III. BRANE INFLATION 



B. Power law potentials 



A. Exponential potentials. 

The simplest case to investigate involves the dynamics 
associated with a pure exponential potential [5], 



V{4>) ~ Vn exp(aiK;</ 
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For completeness we recall the key results. During slow- 
roll inflation, from Eq. (7) we have 
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Inflation ends when e = 1, giving 
Vend — 2a 2 Ah , 
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which implies that typically the term quadratic in the 
density still dominates the linear term at the end of in- 
flation. 

The potential Vn , N e-foldings from the end of infla- 
tion can be obtained from Eq. (10) and is given by the 
simple formula 



V N = Vend (N+l). 
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The amplitude of density perturbations generated dur- 
ing the inflationary period fixes the brane tension 
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which significantly is independent of the mass parameter 

v . 

From Eqs. (3) and (27), we have, H N = (N + l)H end ; 
and substituting into Eq. (12) we find that N = 70.7 for 
A = 0.35, g p = 100 and N = 72.1 for A = 0, g p = 10. 

Inserting this into Eq. (28), and using the observed 
value from COBE of A s = 2 x 10~ 5 [14], we find that 
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indicating a brane tension of order the Grand Unified 
Theories scale. 

For the scalar and tensor spectral indices and their 
derivatives we find, 
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For power law potentials, 

V{4>) = V {Kct>) n 
during slow-roll inflation we have 
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For the particular case of n = —2 the condition for in- 
flation is Ab < Vo/8, i.e. when this condition is satisfied, 
inflation ends in the low energy limit. When n ^ 2, wc 
obtain, at the end of inflation, 
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and N e-foldings prior to the end, 
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The brane tension is no longer independent of the mass 
scale Vq but is related to it by, 
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when n^4, otherwise Ab is completely free and only V 
is fixed, yielding 



V 



75tt 2 



1 



256 k 4 (1 + N) 



-A 2 



(37) 



For the scalar and tensor spectral indices and their 
derivatives we find, 
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which reduce to Eqs. (30) in the limit n — > oo. 

For negative powers of n, the precise number of e- 
foldings corresponding to the scale entering the Hubble 
radius today is easily determined from Eqs. (12) and (35). 

In Fig. 1 we show the dependence of the number of 
e-foldings N with n for four different combinations of 
the time delay between the production of particles and 
the end of inflation and the number of fields experiencing 
particle production. 
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FIG. 1. Dependence of the number of e-foldings on n. 
From top to bottom, for A = 0, g p — 10; A = 0, g p = 100; 
A = 0.35, g p = 10; A = 0.35, g p = 100. 



For completeness we review the main features of power 
law potentials in the case of standard inflation. In this 
case, there is no freedom for the mass scale Vq which is 
fixed by the density perturbations as 



Vn = 
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C. Future observational constraints 

In Fig. 2 we show how a selection of different ob- 
servational parameters depend on the precise value of 
n comparing standard inflation with brane world infla- 
tion. In particular, we consider here the spectral in- 
dex ns, the ratio of tensor to scalar perturbations r = 
AnA^/A 2 ^ = — 27m t and the derivative of the spectral in- 
dex dns/d\nk. For the positive power law case, reheat- 
ing can happen through conventional mechanisms, and 
hence we take the number of e-foldings corresponding to 
the scale entering the Hubble radius today of TV = 50. 
The filled areas correspond to estimated error bars on 
these observables based on upcoming Planck observa- 
tions, according to Ref. [15]. We quote, Sng — 0.006, 
Sr = 0.04 and Sdns/dlnk = 0.01. The error bars are 
centred on a (f> 2 standard inflationary cosmology for the 
sake of illustration. There are several pieces of informa- 
tion one can extract from this figure. Firstly, a <ft 2 theory 
in brane worlds is degenerate with the same theory in 



standard inflation. The case becomes worse for a (j> 4 the- 
ory as the exponential potential presents similar values 
for all the observables. Secondly, in brane worlds, an 
exponential and power law with n < — 10 are degener- 
ate potentials. Thirdly, if we could rule out the expo- 
nential potential, the remaining models would be truly 
non-degenerate only if |n| > 20. And finally, it appears 
that the derivative of the spectral index will be unable to 
give us any information about the underlying inflationary 
cosmology. In summary, based on this small sample, it 
appears that upcoming satellite experiments wont be able 
to conclusively discriminate between the various models 
of inflation. 
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FIG. 2. Dependence of the spectral index, the ratio of ten- 
sor to scalar perturbations and the derivative of the spectral 
index on the parameter n. The solid line: exponential poten- 
tial; dashed line: power law potential in brane worlds; dash 
dotted line: power law in standard inflation. In the case of 
the exponential potential and power law with negative n we 
used Eqs. (12), (27) and (35) with g p = 100 and A = 0.35, 
to define TV. For the power law cases with n positive we used 
TV = 50. 

In the next section, we go one step further trying to 
use the same field to link the early inflation epoch with 
the current accelerated expansion of the universe, which 
is strongly suggested by supernovae type la (SN la) ob- 
servations. 



4 



IV. QUINTESSENCE 



nucleosynthesis and at the time of recombination [19]: 



In the case of reheating by means of gravitational pro- 
duction of relativistic particles, the scalar field does not 
totally decay. Many authors [7,8] have suggested that 
the field could then dominate the energy density at late 
times, providing the dark energy that appears to be driv- 
ing the observed acceleration of the universe today [16]. 

In the brane world picture the potential docs not have 
to have a shallow slope to satisfy the slow-roll conditions 
followed by a steep fall where the latter break, as in stan- 
dard inflation. In fact, we have seen that the potential 
can be exponentially steep and still inflate at early times 
and at the point where the brane correction becomes neg- 
ligible inflation comes gracefully to an end. This fea- 
ture broadens the range of potentials one can consider 
for "quintessential inflation" provided the brane effects 
become negligible before nucleosynthesis [9-11]. 

Using the modified Fricdmann equation (3), Huey and 
Lidsey [9] showed that for a pure inverse power law po- 
tential, the requirement that the field was both sub- 
dominant prior to nucleosynthesis, and tracked the mat- 
ter and radiation in the universe before today implied 
that the value of the power n in Eq. (32), is so large 
(and negative) that the corresponding equation of state 
(w,p = Pcf,/ p<f,) is too large to give an accelerating universe 
at this epoch. They concluded that within the class of 
inverse power law potentials, the only successful model 
would correspond to a kind of "creeping quintessence" , 
in which the energy density of the field is kinetic en- 
ergy dominated from its value at the end of inflation to 
its present value. The field effectively remains frozen 
(with equation of state = —1) until the energy den- 
sity of matter matches its energy density, which occurs 
only around a redshift zwl, At this point the field starts 
rolling down its potential again, and dominates the to- 
tal energy density slowly raising the equation of state 
from — 1, thereby explaining the present day acceleration 
of the universe. Although an interesting scenario, its 
similarity to a pure cosmological constant at late times, 
makes the evolution of this quintessence field very diffi- 
cult to confirm on the basis of upcoming SN la or CMB 
observations, especially if the equation of state of dark 
energy is very close to —1 as observations suggest [17]. 
However, if the contribution of the quintessence field to 
the total energy density is significant, it can increase the 
expansion rate resulting in a shift in the positions and 
an increase in the amplitude of the peaks of the angular 
power spectrum [18], hence opening up the possibility 
that we may be able to discriminate a quintessence field 
from a bare cosmological constant type behaviour. 

In what follows we analyse a class of scalar potentials in 
which the energy density of the quintessence field tracks 
the background energy density during the epoch of de- 
coupling. We are able to constrain the parameters of 
the potentials by employing the values allowed on the 
fractional density of a scalar field, both at the epoch of 



1K UC < 0.045 , 



1K 0C < 0.39 . 



(45) 



In the next subsections we will consider the two possi- 
ble cases of evolution of the scalar field between the end 
of inflation and the epoch the attractor is reached. 



A. Monotonic evolution 

After the initial period of inflation and following that 
of the kinetic energy domination, the scalar field is known 
to freeze at a value given by [20] 
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Using Eq. (6) and the results of the previous section 
for an exponential scalar potential, we can calculate the 
bounds on the allowed slope of the potential. We require 
the field to enter the scaling regime (whereupon fi^ = 
4/ a 2 when radiation is dominant) after nucleosynthesis 
but before radiation-matter equality, i.e., 



v q /" 2 s n^c) < 4 P nuc /« 2 
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These requirements impose the following bounds: 3.6 < 
a < 4.8 for g p = 100, A = 0.35 and 3.4 < a < 4.5 for 
g p — 10, A — 0, results found to be in good agreement 
with our numerical integrations (see row corresponding 
to zero oscillations in Table I). Now, in the standard 
Fricdmann cosmology, a pure exponential potential with 
a slope within this range of parameters has a stable at- 
tractor solution for the energy density that scales with 
the background [21]. Therefore, this solution does not 
provide a late time accelerating universe. A number of 
modifications to this type of potential have been pro- 
posed and investigated in depth which do fit the obser- 
vational data [22-24] . These models present a pure expo- 
nential term dominating for almost the entire history of 
the universe and only recently the other features of the 
potentials become dominant to provide the accelerating 
universe, hence, the above constraints also apply to the 
exponential part of these potentials. In Fig. 3, we show 
the evolution of the equation of state of the scalar field 
for the two exponentials potential (2EXP) [23,10] 



V{(f>) = V (e a ^ + e 13 ^) 



(48) 



We would like to make two remarks on the pure expo- 
nential case. The first concerns the small range of the 
parameter a where this scenario can be implemented re- 
vealing some degree of fine tuning. The second, refers to 
the smallness of the parameter itself. In this range of pa- 
rameters, the contribution of the field is large enough for 
the integrated Sachs-Wolf (ISW) effect to increase the 
power on large scales. Under the assumption that the 
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FIG. 3. The evolution of the equation of state of the scalar 
field after inflation. The dashed line represents the 2EXP 
model with a = 4, j3 = 0.1 and the solid line the SUGRA 
model with n — —23. 

anisotropies we see today are primordial, the normalisa- 
tion of the power spectrum to COBE at the multipole 
£ = 10 therefore, results in a wrong suppression of the 
small scales anisotropies. Hence, for these type of poten- 
tials, the analysis of the angular power spectrum must 
be made with special care to take into account the ISW 
effect. 



B. Oscillatory evolution 

In the case of very steep potentials, the kinetic energy 
of the field at the end of inflation can be large enough to 
overtake the potential by a huge factor. Consequently, 
the field might freeze at a value corresponding to an en- 
ergy density lower than the energy density of the back- 
ground at equality, or even lower than today's energy 
density. One then requires an additional mechanism to 
slow down the field's evolution. Potentials with a mini- 
mum can provide the necessary solution. 

For potentials with a minimum, the field can perform 
several wide damped oscillations around its minimum 
while energy is continually being converted between ki- 
netic and potential. Hence, extra time for the friction 
term in Eq. (4) to slow down the field. In the regime 
of wide oscillations the equation of state of the field is 
close to unity as the potential is very steep and looks like 
V oc 4> n for u» 1. Hence, in practice, the field energy 
density continues decaying as a~ e . In Fig. 3 we show 
the evolution of the equation of state of the field for a 
potential with a minimum (solid line). Fig. 4 shows the 
evolution of the energy densities of the scalar field and 
of the background fluid and the evolution of the scalar 
potential for the same potential with a minimum. After 
inflation the field quickly becomes kinetic energy dom- 
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FIG. 4. The evolution of the energy densities of the scalar 
field (solid) and of the background fluid (dashed) and the evo- 
lution of the scalar potential (dash-dotted) for the SUGRA 
model with n — —23. 

inated. The sharp peaks with w^, = —1 in Fig. 3 or 
the epochs when V(<p) = Kninimum in Fig. 4 represent 
the points at which the field passes through the mini- 
mum. In this example, the field becomes sub-dominant 
just before nucleosynthesis as drops below 0.5 before 
z > 10 10 . After a few oscillations the field has slowed 
down enough and freezes (for z f=a 10 5 , in this example). 
From here the usual dynamics resume. When the back- 
ground fluid's energy density drops, the field reaches the 
attractor solution. The contribution of the field increases 
when it approaches the minimum of the potential, then 
the equation of state quickly decreases to — 1 accelerating 
the universe by today. 

We give two examples corresponding to models pre- 
sented before in the literature. The 2EXP potential with 
a/ (3 < 0, satisfies the above requirements. For the par- 
ticular case of symmetric slopes i.e. [3 — —a [25], we find 
that to satisfy the nucleosynthesis bound and for the field 
to scale after this period and no later than radiation- 
matter equality, the slopes must satisfy the constraints 
presented in Table I. 
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TABLE I. Constraints on tracking quintessential inflation 
for 2EXP with symmetric slopes. 1~ refers to the case in 
which the field freezes when rolling up the potential and 1 + 
to the case in which the field is rolling down before it stops. 



The second example is for the SUGRA type potentials 
[26], 



G 



V(</>) = V (/c^) n exp 



(49) 



For initial conditions of (j> <C ^minimum, the bounds in 
Eq. (45) above place a restriction of — 78 < n < — 23 for 
g p = 100, A = 0.35 and n < -86 for g p = 10, A = 0, 
corresponding to two oscillations of the field (see Fig. 3) 
and w,f, w —0.8, today. 

To obtain the quoted results, the initial conditions of 
the scalar field in our simulations were such that they 
gave a number of e-foldings of inflation consistent with 
the number of e-foldings corresponding to scales entering 
the Hubble radius today in Eq. (12). We have taken 
flc/, w 0.7, today. 



V. CONCLUSIONS 

We have demonstrated, giving examples, of a number 
of classes of scalar potentials that provide realistic models 
of "tracking quintessential inflation" . We have seen that 
the parameter space of potentials described by a pure ex- 
ponential at early times is highly contrived on the basis of 
nucleosynthesis bounds and under the requirement that 
the scalar field had reached the attractor by radiation- 
matter equality. However, potentials with a minimum 
offer a larger range of parameter values for which these 
requirements are satisfied. 

For the potentials we have considered, successful re- 
alisation of "quintessential inflation" requires potentials 
with very steep slopes in order to satisfy the nucleosyn- 
thesis bound, therefore we expect for all of these models 
similar values for the observables, i.e. n$ w 0.94 and 
r w 0.27 (using Eqs. (38) for large n with N rj 70 or 
from Fig. 2). 

We have seen in Sec. Ill and Fig. 2 that these models 
can be degenerate to a <f) 4 theory either in standard or in 
brane world inflation. It is important to know then, if the 
type of potentials considered here fit both the reconstruc- 
tion of the quintessence potential from SN la observations 
[27] and the reconstruction of the inflaton potential from 
CMB observations [28] giving indications of a unified the- 
ory of "quintessential inflation" . 
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